Sub- monolayers of carbon on a— iron facets: an ab initio study 
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Motivated by recent in situ studies of carbon nanotube growth from large transition-metal 
nanoparticles, we study various a— iron (ferrite) facets at different carbon concentrations using 
ab initio methods. The studied (110), (100) and (111) facets show qualitatively different behaviour 
when carbon concentration changes. In particular, adsorbed carbon atoms repel each other on the 
(110) facet, resulting in carbon dimer and graphitic material formation. Carbon on the (100) facet 
forms stable structures at concentrations of about 0.5 monolayer and at 1.0 monolayer this facet 
becomes unstable due to a frustration of the top layer iron atoms. The stability of the (111) facet 
is weakly affected by the amount of adsorbed carbon and its stability increases further with respect 
to the (100) facet with increasing carbon concentration. The exchange of carbon atoms between 
the surface and sub-surface regions on the (111) facet is easier than on the other facets and the 
formation of carbon dimers is exothermic. These findings are in accordance with a recent in situ ex- 
perimental study where the existence of graphene decorated (111) facets is related to increased 
carbon concentration. 

PACS numbers: 31.15.ae,34.50.Lf,36.40. Jn,75.50.Bb,75.70.Rf 



I. INTRODUCTION 

Carbon nanotubes (CNTs) are a versatile material 
with a wide range of potential technological applications 
in fields such as mechanical engineering, electronics and 
biotechnology. The chemical vapor deposition (CVD) 
method has established itself as the most effective way to 
produce CNTs in mass quantities. In this method, car- 
bon containing molecules (hydrocarbons, CO) are disso- 
ciated on catalytic nanoparticles (consisting typically of 
transition metals and their alloys) where carbon eventu- 
ally forms graphitic structures and nanotubes. For wider 
technological exploitation of CNTs, better control over 
the growth process is needed. Controlling nanotube chi- 
rality and preventing particle poisoning and CNT growth 
termination are the most important aspects. In order to 
gain such a control, insight at microscopic level into the 
CNT growth process is needed. This kind of insight can 
be obtained by performing in situ experiments where the 
CNT growth is directly observed. 

Recently, several in situ environmental transmission- 
electron microscopy (TEM) studies of carbon fiber and 
nanotube growth have been carried outP"^. In these 
studies, the catalyst particles were either crystalline 
and/or "liquid-like" (i.e. crystalline with high self- 
diffusivity) during the growth pro cess. S tudied nanopar- 
ticle materials ranged from nickeP^^, cobaltF^^ and 
irorfen^ ^ Q alloys^ of these metals. A common factor in 
many of these investigations is the appearance of step 
edges and new facets as carbon is introduced to the 
nanoparticle and the growth of graphene layers from 
these special regions^EJEJ i n the case of nickel, this phe- 
nomenon was attributed to the stabilization of nanopar- 
ticle step edges upon carbon adsorption and the trans- 
port of catalyst metal atoms away from the step edge 



regionPSEU The energy barrier for carbon bulk diffusion 
in nickel was concluded to be very high when compared 
to any surface related diffusion phenomena.^ 

In several studies, the dominant role of surface and 
sub-surface has been emphasizecP^ while other stud- 
ies, considering mainly iron, suggest the importante of 
bulk diffusion ^ 9 * 12 !. Very recently, CNT growth from car- 
bidic phase (cementite) in iron nanoparticles has been 
demonstrated^^. 

In a very recent study, solid-state a— iron nanopar- 
ticle was encapsulated inside multi-walled carbon nan- 
otubes (MWCNTs) while carbon was injected into the 
nanoparticle by electromigration 6 . The nanoparticle was 
observed to stay solid and crystalline during the growth. 
Similar to earlier in situ studies, new facets, showing 
growth of graphitic material, appeared on the nanoparti- 
cle surface. The orientation of the nanoparticle was ana- 
lyzed and the MWCNT walls encapsulating the nanopar- 
ticle were parallel to the (110) facet. The existence of 
the (111) facet was observed to depend strongly on car- 
bon concentration, and nanotube cap was formed on a 
rounded (100) facet. Ab-initio simulations were per- 
formed for "graphenated" and carbon saturated surfaces 
in order to reproduce the nanoparticle shaped 

In the case of a— iron, the morphology of (110), (100) 
and (111) facets is quite different. This can result in 
very different diffusion barriers, carbon-carbon bond for- 
mation energetics and kinetics. For better understand- 
ing of experiments it is important to perform ab ini- 
tio simulations and correlate computational results to the 
phenomena observed in the in situ studies. In particu- 
lar, surface energies can be used to produce the physical 
shape of the nanoparticle which can be compared with 
the experiments^^. Activation energies of diffusion on 
and into the facets may provide information about rate- 
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limiting steps of graphitic material format iorP^. 

Recently, there have been several ab initio studies re- 
lated to these issues: pure a— iron facets have been stud- 
ied extensively by Blonski and KiejngP^^, while carbon 
adsorption and diffusion on and into the (110) and (100) 
facets were studied by Jiang and Cartel. The (100) 
facet has drawn some attention very recently, as segre- 
gated carbon atoms form stable, periodic structures on 
this facetP^H] Some carbide surfaces have been studied 
with ab initio methods by Chiou and Carter^. 

In this work we study the a— iron (110), (100) and 
(111) facets at different carbon concentrations. We ad- 
dress such topics as the interaction of adsorbate atoms 
in coadsorption configurations, carbon diffusion (for the 
(111) facet) and formation of stable carbon-rich struc- 
tures (mainly on the (100) facet) in the topmost iron 
layer. Relative surface energies as function of carbon 
concentration and energetics of the smallest units in- 
volved in graphitic growth, the C2 molecules, are studied. 
This work is organized as follows: in Sec.(|ll|), simulation 
of iron-carbon systems, different carbon chemical poten- 
tials, calculation of surface e nerg ies and computational 
details are discussed. In Sec. (Ill) the morphology of the 



studied a— iron facets is discussed followed by the com- 
putational results. Discussions and conclusions are made 
in Sec.dlV). 



II. METHODS 



A. Iron-carbon systems 



Iron with dissolved carbon exhibits a complex phase 
diagram as a function of temperature and carbon con- 
centration, where a— iron ("ferrite", bcc crystalline, fer- 
romagnetic), 7-iron ( "austenite" , fee crystalline, antifer- 
romagnetic) and cementite (FesC) are among the com- 
peting phases^!. When considering nanoparticles instead 
of bulk, the phase diagram will be modified; in particu- 
lar, it is known that small iron nanoparticles prefer the 
7-iron phase instead of the bulk a— iron phase^^H. 

While collinear spin calculations within the Density 
Functional Theory (DFT) and General Gradient Approx- 
imation (GGA/ 2 ^ seem to work very well for a— irorP^, 
it is not obvious how to perform calculations for 7-iron, 
as the FCC iron, observed only at high temperatures, 
is paramagnetic. These and some other problems have 
been discussed in detail by Jiang and Carter in Ref.[26 . 

Considering systems involving nanoparticles, it has 
been shown that "large" iron nanoparticles (more than « 
100 atoms) prefer the bcc structure^. This structure was 
also observed in the large iron nanoparticle investigated 
in the recent in situ experimental study^. 



B. Technical Details 

In order to visualize the different a— iron facets, a 
small volume of bulk, cleaved into several directions 
is illustrated in Fig[TJg). The form of this bulk vol- 
ume mimicks the elongated nanoparticles seen in the in 
situ experiment 6 . A real nanoparticle has, of course, sev- 
eral other facets which are not visualized in FigfTJg) and 
not considered in this work. Top and side views of the 
facets are illustrated in panels (a-f) of Figjl] 

We employ the periodic supercell method to model the 
iron facets. An infinite surface is modelled by a slab con- 
sisting of few layers of iron atoms, with sufficient vacuum 
of 15A between the slabs. The slabs we have used are de- 
picted in Fig [I] and they consist of 3x3x6 (the last num- 
ber denoting the number of layers) and 2x2x12 slabs for 
(110), (100) and (111) surfaces, respectively. Areas of 
the unit cells depicted in Figjl] are 51 A 2 , 72A 2 , and 56 
A 2 for (110), (100) and (111) slabs, respectively. 

Next, in order to calculate adsorption and reaction en- 
ergies on the surface, we define a convenient energy quan- 
tity E s (shifted energies) as follows: 



E S (X*) = E(X*) - E 



(1) 



where E(X*) is the energy of the adsorbed surface species 
X* and Eo is the energy of a surface unit cell without 
adsorbates. Now the adsorption energy can be written 
as follows: 



E a 



E(X*) - E(X) -E = E S (X*) - E(X), (2) 



where E(X) is the energy of an isolated atom in vacuum. 

We also present adsorption energies, especially in the 
case of high carbon coverage, by using some other chem- 
ical potential: 

E c = E(X*) - (E + nn c ) = E S (X*) - n/i c , (3) 

where n is the number of carbon atoms and \i c is the 
chemical potential. We take the chemical potential as 
the energy per atom in graphene. Energy E c then re- 
flects the energy cost to accommodate carbon atoms into 
the metal adsorbant instead of graphene. In previous 
wor k^6i20i26j var i ous chemical potentials have been used, 
including the energy per isolated carbon atom and the 
energy per carbon atom in graphite or in graphene. 

The shifted energy values of Eq. can be used to cal- 
culate reaction energetics on the adsorbate, i.e. to look 
at energetics of processes like C*+C* — ^CJ. In the case 
of adsorbate-adsorbate repulsion on the lattice, this can 
provide information about stress release upon dimer for- 
mation on the iron surface. The energy for a reaction 
X*+Y*^XY* can be calculated as follows: 

AE= (E(XY*) + E ) - (E(X*)+E(Y*)), (4) 

This equation can be written, using the energy values E s 
of Eq.Q as follows: 

AE = E S (XY*) - (E 8 (X*) + EsQT*)). (5) 
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Figure 1: (color on-line) Computational unit cells for different a— iron surfaces used in the simulations viewed from (a,c,e) top 
and from (b,d,f) side, (a-b): (110), (c-d): (100) and (e-f): (111) surfaces. Atoms in topmost (lowermost) layers are marked 
with brighter (darker) shades. Atoms in the unit cells used in this work are marked with blue color, (g) A portion of bulk 
a-iron cut at different angles, demonstrating the positions of different crystallographic surfaces. For (100) surface, the topmost 
layer is marked with black colour. For (111) surface, the two topmost layers are marked with blue colour. For (110) and (100) 
the unit cells correspond to 3x3 periodicity, while for (111) to 2x2 periodicity. 



• = octahedral site 




B=1.48 A 

Figure 2: (color on-line) Local coordination at the bulk octahedral site in a— iron and how it is exposed on the (110), (100) 
and (111) facets. The octahedral site is marked in all insets with green color, (a) The coordination of octahedral site in the 
bulk. Coordination of octahedral sites on the (b) (110), (c) (100) and (d) (111) facets. 



In the Results section, we tabulate values of E s and then 
use these tabulated values to calculate reaction energetics 
X* + Y* -^XY* using Eq. g . 

The surface energy G of a specific nanoparticle facet 
with adsorbates has in some earlier works 6 27 been de- 
fined as follows: 

G = E(X*) - NE bulk - n/i c , (6) 

where TV is the number of the surface atoms in the slab 
used for simulations and E bulk is the energy per atom in 
the bulk. As Eq.([6| reflects the energy cost to create a 
surface from the bulk and adsorbing atoms on the surface, 
it is more accurate to use an equation involving explicitly 
the surface energy E sur f as follows: 

G = E sur f + E c , (7) 

where E c is the adsorption energy ofEq.Q. InEq.0 the 
surface energy E sur f is a well defined quantity, while the 



chemical potential of the carbon adsorbates is sensitive 
to the source of carbon atoms. 

In this work we u se the values calculated by 
Blonski and KjiengP 11 ^ for the (110), (100) and (111) 
facets which have been evaluated as discussed by 
BoettgeM These are E snr/ =140 meV/A for the Fe(110) 
and Fe(100) surfaces and E snr f=160 me V/A for the 
Fe(lll) surfaced 

Our calculations were performed in the framework of 
the density fun ctiona l theory (DFT), as implemented in 
the VASP cod d 29 * 3Q [ All calculations were done using 
project or- augmented waves (PAWs^P and the Perdew- 
Burke-Ernzerhof (PBE) generalized gradient approxima- 
tion (GGA)P. We used the Monkhorst-Pack (MP) 
sampling 32 of the Brillouin zone in calculations involving 
the slab. The sampling used was 7x7x1 in the case of all 
the slabs which corresponds to A^^^O.01 A -2 (area in 
the reciprocal space per sampled k-point). A systematic 
search to find the optimal adsorption sites for C atoms 
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and C2 molecules was performed on the slabs of FigjT] 
along the lines of Ref. [33] . 

Spin polarization was included in all calculations. The 
cutoff energy of the plane wave basis set was always 420 
eV. The mixing scheme in the electronic relaxation was 
the Methfessel-Paxton met ho cP^ of order 1. Conjugate- 
gradient (CG) relaxation of the geometry was performed 
and if needed, the relaxation was continued with a semi- 
Newton scheme. This way we were able to reach a max- 
imum force residual of « 0.01 eV/A. In all calculations 
the special Davidson block iteration scheme was used and 
symmetries of the adsorption geometries were not uti- 
lized. 

As carbon chemical potential, we used either the en- 
ergy of an isolated atom in vacuum or the energy per 
atom in graphene. For calculation of the chemical poten- 
tial, identical parameters to those described earlier in this 
section were used. For an isolated, spin-polarized carbon 
atom calculated in a cubic unit cell with 15 A sides, we 
obtained the total energy of E=-1.28 eV. For graphene, 
and using a k-point sampling of 25x25x1 we obtained 
a lattice constant a=2.468A. This is slightly larger than 
values obtained by LP A, ESI but is identical to a previous 
calculation using GGA 36 . The energy per carbon atom 
in graphene we obtained is E=-9.23 eV. 

Nudged elastic band (NEB) calculations^ were per- 
formed with VASP. Atoms in the topmost layer were al- 
lowed to move freely, and in some cases, the atoms below 
the topmost layer were allowed to move into z-direction 
(normal to surface) only. Thus we were able to avoid the 
(artificial) collective movement of the surface slab atoms 
that sometimes occured during the minimization. 



III. RESULTS 
A. Morphology of a— iron facets 

We can expect from some earlier studies concerning 
carbon soluti on into bulk iron and adsorption on iron 
surfaces ) 17 * 26 ' that carbon prefers sites of maximum co- 
ordination: in bulk iron, it prefers the 6-fold octahe- 
dral site^ and on the (110) and (100) facets, carbon 
moves into sites that offer highest possible number of 
iron neighbours^. Keeping this in mind, we will give in 
this section a qualitative picture of carbon adsorption on 
different facets. This analysis is based on the octahedral 
site of bulk a— iron. 

The local coordination of the octahedral site is illus- 
trated in Figj2|a). When carbon is adsorbed into this 
site, a tetragonal distortion in the bcc lattice takes place 
and distance B is expanded. As the bulk is cleaved along 
a specific direction in order to create a surface, the octa- 
hedra become cleaved in a specific way, exposing octahe- 
dral sites. The way the octahedral sites are exposed in 
different facets, has been illustrated in Figj2jb-d). On the 
(110) surface, the exposed octahedral sites have neigh- 
boring iron atoms at distances A and B. In the case of 



the (100) surface, there are several exposed octahedral 
sites where the neighboring iron atoms are simply at a 
distance A. 

Assuming that carbon tries to maximize its coordina- 
tion with iron on the surface (as discussed above), it will 
always prefer an exposed octahedral site, as this kind of 
site offers maximum coordination within the bcc lattice. 
The displacements of iron atoms surrounding a surface 
exposed octahedral should be very similar to the bulk 
tetragonal distortion (i.e. expansion of (B) and a slight 
contraction of (A)). This distortion must be energetically 
very different on the distinct surfaces. Depending on how 
the distortion of A and B fits the facet morphology, quite 
different adsorbate-adsorbate repulsions can be formed; 
for example, expanding B on the (110) facet (Figj2^e)) 
consists of pushing neighboring iron atoms apart. On 
the (100) facet there are many sites with no need to re- 
arrange the iron atoms as only small contraction of A is 
needed. 



B. Bulk a— iron facets 

For the bulk iron lattice constant we obtained 
a=z. 83 A, which agrees well with an earlier computa- 
tional value^ and the experimental value of 2.87 A [38]. 
For the bulk magnetic moment we obtained M=2.18 
jiB • In earlier works, interlayer relax ations for various 
a— iron facets have been studiecP^^. The most impor- 
tant effect is the inward relaxation of the outermost layer. 
The interlayer relaxations can be sensitive to the slab size 
and to the scheme used (symmetric/non-symmetric slab, 
number of fixed layers jflMlsj i n Tabjl] we compare our 
results with previous ones. In our scheme, the slab is 
non-symmetric as the atoms in the three bottom layers 
are fixed. In Ref. [15 three topmost layers were allowed 
to relax, while in Ref. [16 freestanding slabs were consid- 
ered. As evident from Tabjl] we can see that the type 
of relaxation (either expansion or contraction) is quite 
consistent. Magnitudes of expansion/contraction have a 
few differences of the order of 10% in the case of the 
(111) facet, but on the other hand, in this facet the bulk 
interlayer distances are very small (^0.82 A). 

Values of magnetic moment in different layers are pre- 
sented in Tab El Consistent with earlier results^ we ob- 
serve that (100) has the highest top layer magnetic mo- 
ment and that in all slabs, the value of magnetic mo- 
ment approach to that of bulk as we move inside the 
slab. Our values for the moments in the topmost layer, 
2.56/iB (HO), 2.96/iB (100) and 2.84/i B (111), compare 
well with the values of Ref. [16], namely 2.59/i# (HO), 
2.95/iB (100) and 2.81/i B (111). 

C. Atomic carbon, coadsorption, dimers 

In this section we study atomic carbon adsorption (l/9 
ML coverage for (110) and (100) and 1/4 ML coverage 
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0Ll5 




d 6 7 




(110) 


-0.1 


0.3 


-0.5 


-0.2 0.04 


0.2 


Ref. [16] 


-0.11 


1.16 


1.14 








Ref. [15] 


-0.4 


0.5 


-0.7 








This work 


(100) 


-3.6 


2.3 


0.4 


-0.4 


-0.01 


-0.5 Ref.fE] 


-3.09 2.83 


1.93 








Ref. [15] 


-1.2 


3.4 


3.5 








This work 


(111) 


-17.7 


-8.4 


11.0 


-1.0 


-0.5 


0.1 


Ref. [16] 


-6.74 


- 16.89 


12.4 








Ref.[l5] 


-3.8 


-18.0 


11.0 0.0 


-0.7 


1.6 


This work 



Table I: Interlayer relaxations in the slabs of Fig^as percent- 
age of the bulk distances, dij is the distance between layers i 
and j. 





1 2 3 4 5 6 


110 


2.57 2.28 2.17 2.17 2.29 2.56 


100 


2.96 2.37 2.47 2.50 2.38 2.96 




1 2 3 4 5 6 


111 


2.89 2.39 2.49 2.30 2.28 2.15 


7 8 9 10 11 12 




2.20 2.24 2.29 2.46 2.29 2.84 



Table II: Magnetic moment of atoms in different layers (/jb ) 
in the slabs of Fig[l] Number 1 denotes the bottom layer. 

for (111)) as well as dimer- and co- adsorption. We study 
in detail how the adsorbates either repel or attract neigh- 
boring iron atoms in the topmost iron layers and will use 
the resulting displacements of iron atoms as our leading 
argument when describing the energetics at higher car- 
bon concentrations in the next section. 



Adsorbate|E ads (eV)|E c (eV)|E a (eV)|BL (A) 


Fe(110) 


C-l 


-7.98 


-0.03 


-9.26 




C-2 


-6.91 


1.04 


-8.19 




C-3 


-5.48 


2.47 


-6.76 




C 2 -l 


-8.19 


0.73 


-17.72 


1.35 (1.31) 


C 2 -2 


-8.04 


0.88 


-17.57 


1.38 


C2-3 


-7.12 


1.81 


-16.64 


1.32 


CA-1 




0.2 


-18.26 




CA-2 




0.74 


-17.71 




Fe(100) 


C-l 


-8.45 


-0.5 


-9.73 




C-2 


-7.2 


0.74 


-8.48 




C-3 


-7.18 


0.77 


-8.46 




C2-I 


-8.26 


0.67 


-17.79 


1.33 


C2-2 


-7.91 


1.01 


-17.44 


1.36 


C2-3 


-7.51 


1.42 


-17.03 


1.44 


CA-1 




-1.02 


-19.48 




CA-2 




-0.92 


-19.38 




Fe(lll) 


C-l 


-7.74 


0.2 


-9.02 




C-2 


-7.69 


0.26 


-8.96 




C-3 


-7.43 


0.52 


-8.71 




C-4 


-7.41 


0.53 


-8.69 




C2-I 


-8.95 


-0.02 


-18.47 


1.40 


C2-2 


-8.87 


0.05 


-18.4 


1.38 


CA-1 




0.61 


-17.84 




CA-2 




0.82 


-17.63 





Table III: Adsorption energies E a d s (see Eq.pV), energies E c 
(see Eq.(|3|) and shifted energies E s (see Eq7|l])). Values of 
E s can be used directly to calculate reaction energies on the 
surface by using Eq.|5|. Values for C atoms and C2 molecules 
in different adsorption geometries on the a-iron (110), (100) 
and (111) surfaces have been tabulated. Bond lengths (BL) 
for adsorbates and in vacuum (in parenthesis) are listed. Sites 
and geometries have the same labels as in Figs(3|4] and m 
Tab |IV| Coadsorption geometries are tagged with the label 
"CA". 



1. Atomic carbon 

Adsorption geometries for carbon atoms (C-l, C-2, 
etc.) on the different facets are illustrated in Fi gs [3} 
[5] Corresponding energetics are tabulated in Tab jlll) 
Taking a closer look at Figs(3][5j we can see that carbon 
favors the ex posed octahedral sites of Fig(2] as discussed 
in Sec.dfflA). 



The most favorable (C-l) sites fo r (110) and (100) are 
consistent with previous calculations™^ From Tab |m] 
the adsorption energies are -7.98 eV and -8.45 eV for the 
(110) and (100) surfaces, respectively. These compare 
well with earlier computed values of -7.92 eV[17] and - 
7.97 eV© for (110) and with -8.33 eV[17 , -8.335 eV[20] 
for the (100) surface. 

The iron atom below the adsorbed carbon atom on 
100/C-l shifts downwards, corresponding to expansion 
of (B) in Fig 1 2} The coordination of carbon on (100) 
is fivefolcff^^and it is bonded to the iron atom below 
at a distance of 1.98 A. In the following we analyze in 
more detail the intralayer relaxations in the topmost sur- 



face layer, using as a guide th e qual itative discussion of 
these relaxations made in Sec. (Ill A); this kind of analy- 
sis, based on the tetragonal distortion of bulk adsorption, 
was made to some extent in Ref . [17] , but only for the case 
of sub-surface adsorption. 

In 110/C-l of Fig|3] there are considerable intralayer 
relaxations in the topmost layer. The distance between 
iron atoms (b) and (d) contracts by 5 % (0.20 A, corre- 
sponding to contraction of A in Figj2| while the distance 
between (a) and (c) expands 23% (0.65A, corresponding 
to expansion of B). In 100/C-l of Fig J4] the displacement 
of iron atoms is smaller: now both distances (bd) and 
(be) contract only by 5 % (0.15 A), corresponding simply 
to the slight contraction of A; due to the specific cleav- 
ing of the octahedron in Fig{2j the Fe(100) surface offers 
exposed octahedral sites for carbon with very little need 
to move the surrounding iron atoms. 

From Tab |III| we can see that adsorption into 100/C-l 
is 0.47 eV more favorable than adsorption into 110/C- 
1 (similar to the value of 0.41 eV obtained in Ref. [17]). 



Fe(llO) 






Figure 3: Some of the most stable geometries for C and C2 on the Fe(llO) surface. Different geometries are tagged with the 
same labels as in Tab |III| Coadsorption geometries, where atoms are adsorbed into the same unit cell are tagged with the label 
"CA". Iron atom displacements upon carbon adsorption in (C-l) have been marked with arrows. These can be related to Fig[2] 
In (CA-2) one adsorbed carbon atom is pushed towards vacuum. 



Fe(100) 




Figure 4: Some of the most stable geometries for C and C2 on the Fe(100) surface. Different geometries are tagged with the 
same labels as in Tab |III| Coadsorption geometries, where atoms are adsorbed into the same unit cell are tagged with the label 
"CA". Iron atom displacements upon carbon adsorption in (C-l) have been marked with arrows. These can be related to Fig[2] 



When looking at energies E c , we observe that at a low 
coverage of 1/9 ML (i.e. single adsorbed atom), it is 
more favorable for the carbon atom to be adsorbed into 
the iron surface than to be incorporated into graphene. 
While in the case of (110), this tendency is very weak 
(only ~ 30 meV), for (100) it is more significant (0.5 
eV). As will be discussed below, E c is very sensitive to 
the amount of carbon adsorbed on the facets: at lower 
concentrations than we are considering in this paper (less 
than 1/9 ML), E c should clearly become negative also for 
(110). 

The remaining adsorption geometries for (110), i.e. 



110/C-2 and 110/C-3 are metastable down to w 0.012 
meV/A and they lie more than 1 eV higher in energy than 
110/C-l. Their characteristics (local minimum, higher 
order saddle point, etc.) have been discussed in more 
detail in Ref . [T7] . . 

Comparing Figj2|d) and the optimal adsorption sites 
of carbon atoms in Fig(5]we can see that carbon prefers 
exposed octahedral sites on the (111) facet. In 111/C-l, 
there is a 0.74 A expansion in the distance of atoms (b-d) 
and a 0.2 A contraction in the distance of atoms (a-c), cor- 
responding again to (B) and (A) in (Figj^a)). While ad- 
sorbate lll/C-2 exhibits very similar distortions, 111/C- 
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Figure 5: Some of the most stable geometries for C and C2 on the Fe(lll) surface. Different geometries are tagged with the 
same labels as in Tab |III[ Coadsorption geometries, where atoms are adsorbed into the same unit cell are tagged with the label 
"CA". Iron atom displacements upon carbon adsorption in (C-l) have been marked with arrows. These can be related to Fig[2] 



3 breaks the trend a bit as it does not adsorb into an 
exposed octahedral site. It finds a high coordination by 
moving atoms (a),(b) and (c) instead. Atoms (a),(b) and 
(c) all move symmetrically ~0.3A and their distance to 
the carbon atoms becomes 2.1 A. There is also one iron 
atom directly below the carbon at a distance of 1.85 A. 
Adsorbate lll/C-4 is simply a carbon atom adsorbed at 
a bulk-like octahedral site. 

When looking at the adsorption geometry C-l, we can 
observe that it is by definition a "sub-surface" site, i.e. 
the carbon atom resides below the topmost iron layer. 
On the other hand, it has not yet obtained a coordination 
with surrounding iron atoms similar to that in bulk. On 
the contrary, C-2 is clearly a "surface" adsorption site. 
The energy difference between C-l (a "semi" sub-surface 
site) and C-2 ("surface" site) is minimal, only 50 meV. 

The energetics for carbon adsorption on (111) in 
Tab |III| are not directly comparable to those reported in 
Ref.[6j, as in that work the motion of iron atoms was 
constrained (some of the sites will have diffent local ge- 
ometries upon relaxation). Our values for the C-l (-7.74 
eV) and C-3 (-7.43 eV) sites are very close to the value 
reported in Ref.[6j (-7.60 eV) for a similar site. 



2. C2 dimer 

Assuming that carbon atoms in the dimer prefer sim- 
ilar high-coordinated sites as the individual atoms while 
maintaining a reasonable carbon-carbon bond length, 
there are not good possibilities to achieve this on the 
(110) and (100) surfaces, as the optimal C-l sites lie far 
away from each other. The case of (111) is very different; 
looking at Fig{2|d) we can see that there is an abundance 
of optimal adsorption sites within the bond length of a 



carbon dimer. We can then expect that C2 dimer is most 
stable on the (111) surface. 

The optimal C2 adsorption geometries on the (110) and 
(100) of Figs{3]4]are in both facets quite similar: individ- 
ual carbon atoms are 1-3 fold coordinated to iron: one 
Fe-C distance in both cases is 1.85 A, while the remain- 
ing two Fe-C distances are « 2.0A. The C-C bond length 
for IIO/C2-I is slightly expanded while for IOO/C2-I it is 
closer to the isolated C2 bond length. In IIO/C2-I, both 
iron atom distances ac and bd expand by ~5% while 
in IOO/C2-I, the distance between a and c is expanded 
by ~15%. Of the remaining C2 adsorption geometries 
IIO/C2-2 exhibits similar trend as C2-I: both C atoms 
reside in a site ~3-fold to iron atoms. Other C2 adsorp- 
tion geometries are trying to adopt positions where the 
adsorbed carbon can reside in C-l, C-2 or C-3 sites 

The optimal position, C2-I for the carbon dimer on 
(111) is depicted in FigjE] Both carbon atoms reside in 
a C-l site while the Fe-C distances are ~2 A. The C-C 
length in the dimer is expanded by «6%. The adsorbate 
C2-2 exhibits a similar trend. Comparing the adsorp- 
tion ener gies of C2-I in the case of the different facets 
(Tab |III| ), the adsorption of the carbon dimer on 111 is 
at least « 0.7 eV more favorable than on the other facets; 
as discussed above, this is because the (111) facet offers 
nearby optimal adsorption sites for the carbon atoms. 

The reaction energetic of Tab |IV| further demonstrate 
that it is much more favorable to form carbon dimers on 
the (111) facet than on the (110) and (100) facets. The 
(100) facet favors less dimer formation than the other 
ones, as having carbon in atomic form is energetically 
most favorable on this facet. 
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3. Coadsorption 

We have studied coadsorption configurations by plac- 
ing two carbon atoms in a sublattice of the most optimal 
adsorption sites of the individual carbon atoms. The op- 
timal coadsorption sites we found are marked with tags 
"CA" in Figs |3pl 

In the case of (110) (Figj4| the two coadsorption config- 
urations considered become very different; in 110/CA-l, 
both atoms reside in a 110/C-l adsorption site. However, 
in 110/CA-2, one of the atoms is forced to move from the 
C-l site towards vacuum. As discussed earlier, this re- 
sults from the expansion of distance (B) (Figj2| on the 
(110) facet; in 110/C-l, atoms (a) and (c) are pushed 
apart and in 110/CA-2 the carbon atoms are pushing 
the same atom (c) into opposite directions. This creates 
strong repulsion between the carbon atoms and only one 
carbon atom can be accomodated into the C-l site. This 
results in a notable, 0.5 eV energy difference between the 
110/CA-l and 110/CA-2 configurations 

The situation is very different on the (100) surface. As 
discussed above, the adsorption of carbon to the 100/C-l 
site does not involve considerable motion of the surface 
iron atoms, as the only displacement needed is the very 
small contraction of (A) (Figj2|. In the optimal coadsorp- 
tion geometries 100/CA-l and 100/CA-2 of Fig|] there 
are indeed very minor displacements of iron atoms to- 
wards the adsorbed carbon. In CA-2 both carbon atoms 
attract iron atoms (c) and (d), while in CA-1 they pull 
only one common iron atom (c). In Tab. Ill we can see 



that this results in a small 0.1 eV energy difference be- 
tween CA-1 and CA-2. 

In the case of the (111) surface, there are quite many 
neighboring optimal adsorption sites (C-l and C-2) and 
so the number of coadsorption configurations becomes 
large. The two most optimal configurations we found 
(CA-1 and CA-2) are illustrated in FigjI] As evident 
from Tab|III[ their energy difference is only 0.2 eV. 



D. Higher carbon concentrations 

In this section, we study the effect of high carbon con- 
centrations (>2/9 ML for (110) and (100), >3/4 ML for 
(111)) on energies E c (Eq.([3|). Carbon is adsorbed on 
the sublattice formed by the most favorable C-l adsorp- 
tion sites. Our study is most systematic for the (100) 
surface as there is an obvious way how to place an in- 
creasing number of carbon atoms on the surface: we do 
not expect considerable displacement of iron atoms from 
the bulk positions nor the displacement of carbon atoms 
away from the optimal C-l adsorption sites. On the other 
hand, as we saw in previous sections, on (110) surface ad- 
sorbed carbon atoms start to repel each other, while on 
(111) we can expect dimer formation. Our objective for 
(110) and (111) is simply to further demonstrate these 
points (adsorbate repulsion, dimer formation) at higher 
carbon concentrations. We start with the case we studied 



most systematically, i.e. with the (100) surface. 

(100) surface. Carbon concentration of up to 1 ML 
has been considered for the (100) facet in Figj6] A mini- 
mum of energy E c as function of carbon concentration oc- 
curs at the coverage of 6/9 = 0.667 ML at an adsorption 
configuration "3^/(2) x \/(2)" that has been reported ear- 
lier in an experimental^ and theoretical study; in this 
adsorption configuration, carbon atoms form an infinite, 
two atom wide ribbon on the C-l sublattice. There exists 
even more favorable coadsorption configura tion "c( 2x2)" 
which takes place at the coverage of 0.50 ML l 17 l 18 l 2Q [ Both 
of these stable structures have been analyzed using STM 
and/or LEED^ 11 ^ For the 6/9 ML coadsorption struc- 
ture presented in Fig(6] we observe same phenomena as 
reported in an earlier study 20 , most notably the displace- 
ments of iron atoms in the topmost layer as illustrated 
schematically in Figj6] 

The most interesting phenomena in the case of the 
(100) surface is the behaviour of E c in Figj6] as func- 
tion of carbon concentration: the energy lowers as more 
carbon is adsorbed, implying the existence of stable iron- 
carbon phases, while at increased >0.5 ML carbon con- 
centrations the situation becomes unstable and graphene 
formation is favored. This behaviour can be explained 
by considering the displacements of iron atoms in the 
topmost surface layer. 

From the reconstruction patterns presented by arrows 
in Fig (6] we can see that iron atoms move always towards 
the adsorbed carbon. These changes in iron atom po- 
sitions are typically in the range of ~0.2-0.3A and, this 
is in accordance what was discussed in earlier sections. 
As carbon concentration increases, iron atoms become 
increasingly "frustrated" as they are surrounded by car- 
bon; at 1.0 ML, an iron atom has so many neighboring 
carbon atoms that it cannot obtain optimal bonding con- 
ditions with any one of them. 

This frustration mechanism results in sudden change in 
the energetics on the (100) facet; even at such a high con- 
centration as 8/9 ML, it is more favorable to adsorb car- 
bon atoms to the C-l sublattice than incorporate them 
into graphene, but when going from 8/9 ML to 1 ML, 
this tendency is suddenly reversed. In the last row and 
column of FigjHJ we present the energy gain when releas- 
ing surface frustration by dimer formation. It is quite 
large, over 1 eV. 

(110) surface. Some configurations of carbon for 
up to 1 ML have been considered for the (110) facet in 
Fig(7| The initial configurations consisted of adsorbing 
carbon atoms on either the C-l, C-3 or C2-1 sublattice 
(see Figj3l) and the final, relaxed geometr ies are presented 
in Figj7[As we discussed in Sec. QUI C 3| ), at 2/9 ML cov- 
erage there are coadsorption configurations where carbon 
adsorbates start to repel each other. Such effects be- 
come more important at higher carbon concentrations. 
In Fig {7] and for a 3/9 ML concentration, there is still 
one possibility to adsorb carbon without creating an ex- 
cess of surface stress, similar to CA-1 of Fig(3] and this 
is illustrated in the string-like configuration of C3-1. If 
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Fe(llO) 


Fe(100) 


Reaction 


AE 


Reaction 


AE 


2 (C-l) -+ C 2 -l 
2 (C-2) C2-I 
(C-l) + (C-2) C2-I 


0.79 
-1.35 
-0.28 


2 (C-l) -> C2-I 
2 (C-2) C2-I 
(C-l) + (C-2) C2-I 


1.67 
-0.82 
0.43 


CA-1 C2-I 
CA-2 C2-I 


0.53 
-0.01 


CA-1 C2-I 
CA-2 -> C2-I 


1.69 
1.59 


2 (C-l) -> CA-1 
2 (C-l) CA-2 


0.26 
0.81 


2 (C-l) -> CA-1 
2 (C-l) CA-2 


-0.02 
0.08 



Fe(lll) 


Reaction 


AE 


2 (C-l) -> C2-I 
2 (C-2) C2-I 
(C-l) + (C-2) C2-I 


-0.42 
-0.54 
-0.48 


CA-1 -> C2-I 
CA-2 C2-I 


-0.63 
-0.84 


2 (C-l) CA-1 
2 (C-l) CA-2 


0.21 
0.41 



Table IV: Reaction energies for forming C2 dimers and for b ringi ng carbo n ato ms into coadsorption configurations. Geometries 
are tagged with the same labels (C-l, C-2, etc.) as in Tab[IIl]and Figs(3|4] Reaction energies are calculated by taking the 
corresponding energies E s from Tab III and applying Eq.|5]). 
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Figure 6: Adsorption geometries where carbon atoms are adsorbed in the sublattice formed by the C-l sites (see FigQ. Nc 
is the number of carbon atoms and N con / the number of all possible coadsorption configurations when Nc atoms are adsorbed 
on a (3x3) C-l sublattice. For a particular coadsorption geometry illustrated in the table, energy E c and the mean Bader 
electron occupations for carbon atoms (Qb) have been tabulated. Displacements of some iron atoms upon carbon adsorption 
have been highlighted by blue arrows. In the last row and column, the energy gain when forming a C2 molecule at the 1 ML 
concentration is calculated. 



we place atoms at the same 3/9 ML concentration to a 
C3-1 sublattice, only one carbon atom is adsorbed into 
a C-l site and the remaining atoms form spontaneously 
a dimer. This is demonstrated in configuration C3-2 of 
FigJ?) 

At 4/9 ML coverage we have considered several differ- 
ent configurations: C4-1 is very similar to CA-2 of Figj4] 
and the surface stress is again released by expulsing two 
atoms towards vacuum while maintaining at least two 



other atoms on the C-l site; in C4-4, a complex recon- 
struction, where all carbon atoms are residing in 3-fold 
adsorption sites occurs; at 6/9 ML concentration, even 
iron atoms are expelled from the surface in order to re- 
lease surface strain and to accommodate carbon atoms in 
the C-l site; in C6-2, a C-Fe-C molecule is spontaneously 
formed. 

(Ill) surface. In the case of the (111) we discussed 
how there is an abundance of open octahedral sites (C-l 
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Figure 7: Relaxed adsorption geometries on (110) after carbon atoms have been placed in a sublattice formed by the C-l 
(C3-1, C4-l,C6-2), C-3 (C3-2, C4-4, C6-1) and C 2 -l (C4-2, C4-3) configurations of Fig(3] Number of carbon atoms per unit 
cell (Nc) and energies E c have been tabulated. 
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Figure 8: Relaxed adsorption geometries on (111) after carbon atoms have been placed in a sublattice formed by the C-l and 
C-2 sites (see Fig|5}. Number of carbon atoms per unit cell (Nc) and energies E c have been tabulated. 
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Figure 9: Surface energy G (Eq.|7])) as function of carbon coverage (atoms per A 2 ) on the (110), (100) and (111) facets. 
Energies correspond to the fully relaxed adsorption geometries of Figs |3|8| Energy at zero concentration has been fixed to the 
surface energy of the facet (l40 meV/A 2 for (110) and (100), 160 meV/A 2 for (111)), as described by Eq-jT]). 



and C-2 in Figj5|. We also noted that these sites are close 
to each other so that C2 dimers are naturally formed. 
Some of the most favorable coadsorption configurations 
at higher concentrations on the (111) surface have been 
illustrated in FigjHJ From the energetics in FigjHJ we see 
that there is an energy cost when adsorbing more carbon 
on the surface, but not as high as in the case of (110); 
no strong repulsion of carbon atoms takes place as there 
are plenty of adsorption sites and the top layer atoms 
in (111) are flexible to move. In (110), the C2 dimers 
were formed above the topmost surface, resulting from 
expulsion of carbon atoms from the C-l sites, while in the 
(111), the carbon dimers are formed below the topmost 
layer, in the very same sites where the atomic carbon is 
adsorbed. 



Surface energies. The energetics of Tab |III| and 
Figs |6p are also shown in Figj9]as surface energies (em- 
ploying Eq.Q). This figure is not to be taken as an 
exact representation of surface energies as function of 
carbon concentration; situations with "graphenated" sur- 
faces (see Ref.[6]) and many possible carbon adsorption 
configurations are missing. However, one can see some 
clear trends: the steep rise in the surface energy as func- 
tion of carbon concentration in the (110) facet implies 
aggressive graphene formation on this facet; the (100) 
facet forms stable carbidic phases near 0.5 ML = 0.08 
1/A 2 ; at high carbon concentrations the (111) facet gains 
in relative stability with respect to (100) facet. This sta- 
bilization can be understood in terms of frustration of 
the (100) facet, and from the bigger number of adsorp- 



tion sites available and the flexibility of the topmost iron 
atoms on the (111) facet, as discussed above. 



E. Diffusion 

Diffusion on and into the (110) and (100) facets has 
been described earlier by Jiang and Cartel 17 . We re- 
peated their calculations for the activation barriers of 
carbon diffusion on the topmost iron layer and obtained 
identical results. Here we report resuls for carbon diffu- 
sion on the (111) facet which, to our knowledge, have not 
been calculated earlier using ab initio methods. 

When studying sub-surface adsorption and diffusion, 
one must keep in mind the deep interlayer relaxations 
of pure iron slabs and the elongation/contraction of dis- 
tances (A) and (B) (Figj2| which might propagate far 
in the lattice. Then a thorough investigation of activa- 
tion energies and energy cost for carbon atom to enter 
the bulk would require calculations with very thick slabs, 
considering increasing adsorption depths until bulk val- 
ues are recovered. 

In Ref.[17 barriers and energetics for (110) and (100) 
were calculated down to the first sub-surface layer. As 
discussed in Sec.(IIICl) for (111) facet, the 111/C-l site 
of Fig(5] can be classified as "semi" sub-surface site. In 
our calculations, we have considered also the first "true" 
sub-surface site (i.e. a site that has bulk-like coordina- 
tion), but have not pursued the calculation of the diffu- 
sion barrier when going very deep inside the slab as the 
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1.12 eV 





Figure 10: (a) Diffusion on (111) into the C-l and C-2 sites (FigQ and (b) into a deeper sub-surface site. 



limits of our finite slab are quickly reached. 

Carbon diffusion between neighboring C-2 sites and 
between the C-l and C-2 sites has been illustrated in 



Fig 10 'a). Diffusion from C-2 to deeper inside the slab 
and reaching a site where carbon has similar coordination 
as in the bulk octahedral site, is illustrated in FigflO^b). 
This site, not reported in FigJFJ should not be confused 
with the C-4 adsorbate which lies even deeper inside the 
slab. 

The activation energy E a =1.12 eV for surface diffusion 
reported here for the (111) surface is smaller than for the 
(100) surface (1.45 eV), but slightly bigger than for the 
(110) surface (0.96 eV)[T7]. Going from C-2 site to the 
"semi" sub-surface site C-l has a very low barrier of 0.43 
eV. Going deeper inside the slab has a barrier of 0.77 eV. 
As illustrated in the insets of Fig 10 'b), in the transition 
state of the minimum energy path, one atom in the top- 
most iron layer is pushed away from the carbon atom. 
However, the barrier is low (0.77 eV) because the iron 
atoms in the (111) topmost layer are low coordinated 
and flexible to move. We identify the 0.77 eV barrier 
tentatively as the activation energy for carbon atom to 
enter the surface and it is smaller than those reported 
for (110) (1.18 eV) and (100) (1.47 eV). Moreover, from 
Tab |III| the energy difference between the 111/C-l and 
lll/C-4 sites, is only 0.33 eV, which is smaller than the 
same energy difference in (110) (0.62 eV) and (100) (1.19 
eV)[T7]. There are then clear indications that the ex- 
change of carbon atoms between the topmost and deeper 
layers is easier on the (111) facet than on (110) and (100) 
facets. 



IV. DISCUSSIONS AND CONCLUSIONS 

This work was motivated by the recent in situ studies 
of carbon nanotube growth from "large" iron nanopar- 
ticles where different nanoparticle facets seem to behave 



in a very different manner. The facets studied in this 
work correspond to those identified in Ref. [6 , namely the 
a— iron (110), (100) and (111) surfaces. We have studied 
the effect of adsorbing increasing amounts of carbon on 
these surfaces 

The repulsion among adsorbed carbon atoms on the 
(110) facet is strong and because of this, carbon atoms 
are expulsed from the optimal adsorption sites when car- 
bon concentration is increased. This may happen already 
at the relatively low 0.22 ML concentration, resulting in 
dimer and graphene formation. 

The (100) facet behaves in a very different way, as 
stable, carbon rich structures are formed near 0.5 ML 
concentration. When approaching 1 ML carbon concen- 
tration, (100) is destabilized due to frustration of the top 
layer iron atoms, which can be released by dimer and 
graphitic material formation. 

The (111) surface behaves again in a distinct manner; 
the surface energy is rather insensitive to the amount of 
carbon adsorbed on it (at similar atoms/ A 2 concentra- 
tions where (100) becomes unstable). This results from 
the abundance of adsorption sites and from the flexibility 
of the topmost iron atom layer. 

Carbon nanotube and graphene formation are known 
to be highly kinetic processes, so one must be cautious 
when relating the present work - based mainly on to- 
tal energies of different adsorbed carbon concentrations 
- to these processes. However, the arguments concerning 
the nature of adsorbate repulsions and surface iron atom 
frustration on the different facets are valid also in a dy- 
namical situation and at higher temperatures, as long as 
the nanoparticles are crystalline. 

The sudden carbon "supersaturation" at the (100) 
facet near 1 ML carbon concentration might be related 
to the lift of graphitic caps from this facet as more car- 
bon is injected into the nanoparticle as seen in the in 
situ experiment^. When carbon feedstock is exhausted 
in the experiment, the (111) facets start to shrink, even- 
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tually disappearing. This is very likely related to our 
computational results: the relative stabilities of (100) 
and (111) were observed to be sensitive to the amount 
of adsorbed carbon. At low carbon concentration (100) 
gains in stability while at high concentrations, (111) be- 
comes equally stable. The exchange of carbon between 
the topmost iron layer and the subsurface layers was seen 
to be easier on (111) than on the other facets. Carbon 
dimer formation was observed to be most favorable on 
(111). This could be related to the typical TEM image of 
a CVD grown multi-walled carbon nanotubes that show 
several graphitic layers emerging from the (111) faceiP. 

These last observations can be important when trying 
to understand the growth mechanisms of carbon nan- 
otubes in general: in a pure a— iron nanoparticle, the 
portion of (111) facets is very small 16 . As carbon con- 
centration on the nanoparticle surface gets higher, (111) 
facets will be stabilized. Once these facets have been es- 
tablished, graphitic material growth from them can pro- 
ceed as they favor C2 formation and as the movement of 
carbon atoms between sub-surface and surface is easier. 
However, in order to make this idea more solid, more in- 
vestigation about the diffusion of carbon and kinetics of 
graphene growth on the (111) facet must be performed. 



To summarize, carbon concentrations of up to one 
monolayer were studied on a— iron facets. Such aspects 
as repulsion between carbon adsorbates on the (110) the 
frustration of iron atoms on the (100) surface and the 
dimer formation on (111) facet together with their effect 
on the surface energies were discussed. Diffusion on and 
into the (111) facet was studied. Our findings were re- 
lated to a recent in situ study where the appearance 
of (111) facets correlates with increased carbon concen- 
tration. A general idea where increased carbon concen- 
tration stabilizes the (111) facets followed by growth of 
graphitic material from these facets was proposed. 
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